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ABSTRACT 

We calculate the temperature and ionization balance in an outflow from an 
accreting black hole under illumination by hard radiation from the central object. 
Electron scattering of the Fe K a photons within the highly ionized expanding flow 
leads to a decrease of their energy (redshift) which is of first order in v/c, when 
v is the outflow velocity (v is much less than the speed of light c). This photon 
redshift is an intrinsic property of any outflow for which divergence is positive. 
We also find that the equivalent widths of red skewed Fe Ka originated in the 
wind is order of keV. We conclude that redshifted lines are intrinsic properties 
of the powerful outflows that are observed in many compact objects. Downscat- 
tering of the primary line photons generated in the outflow (a more natural and 
probable mechanism than the general relativistic effects in the innermost part of 
the accretion flow) leads to the formation of red-skewed lines. 

Subject headings: accretion — stars: radiation mechanisms: nonthermal — black 
hole physics — atomic processes — line: format ion — galaxies: active 

1. Introduction 

There is a great need for a coherent physical picture predicting how X-ray spectra are 
produced near galactic and extragalactic black hole (BH) sources. Elvis (2003) points out 
that the BH-disk-jet theory is presently a "high" theory, dealing only with overall energetics, 
and describes a naked object devoid of any of the veiling gas that presumably creates the 
observable spectral continuum and atomic features. This "high" theory fails to explain the 
various correlations between timing and spectral features observed in BH sources [for exam- 
ple, see details of these correlations in Mauche (2002); Titarchuk & Wood (2002); Titarchuk 
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& Fiorito (2004)]. Elvis concludes that the community needs a "low" theory with sufficient 
detail to predict the emission (and presumably timing) and absorption phenomenology of 
BHs. 

Many authors [e.g. Tanaka et al. (1995); Nandra et al. (1997); Wilms et al. (2001); 
Uttley et al. (2003); Miller et al. (2004)] have found unusual curvature (red-skewed features) 
at energies between 2 and 8 keV in X-ray spectra of a number of galactic and extragalactic 
BH sources. Significantly, the fluxes of these features remain nearly constant despite the 
large changes in the continuum flux from the central object as shown in the observations 
of MCG-6-30-15, NGC 4051 and others [Markovitz, Edelson, & Vaughan (2003), hereafter 
MEV03], suggesting an origin in the outer regions of the accretion flow. In fact, MEV03 
discuss a possibility of explanation of the lack of correlation between the continuum and 
the line using the ionized disk skin model by Nayakshin (2000) and Nayakshin & Kallman 
(2001) (see also Ross, Fabian & Young 1999). This model is consistent with the prediction 
that both broad and narrow K Q lines should track continuum variations on time scales of 
months to years. However such correlated behavior is not seen through the entire sample. 

In this paper we further discuss a model for redshifting of emission lines by repeated elec- 
tron scattering in a diverging outflow or wind. The basic mechanism is described by Laurent 
& Titarchuk (2004). It is distinct from the proposals for redshifts previously associated with 
outflows, Comptonization (in a low temperature electron cloud, see Sunyaev & Titarchuk 
1980, for details) and the general accelerating outflow model of Sobolev (1957), both of which 
were correctly ruled out by Fabian et al. (1995) in connection with the redshifted Fe Ka 
in MCG-6-30-15. Comptonization would produce an unobserved spectral break at energy 
E ~ m e c 2 /r| ~ 20 keV where r T is the outflow Thomson optical depth, The accelerating 
outflow model required optically thick clouds moving with velocities v ~ 100 000 km s -1 . 
They argue that the major problem with this outflow model is that MCG-6-30-15 has a 
warm absorber. The OVII edge in the absorber indicates that the flow along line of sight 
is less than about 5000 km s" 1 (Fabian et al. 1994), which seriously conflicts with outflow 
of about 100 000 km s -1 . Laurent & Titarchuk (2004) confirm that velocity of the warm 
absorber outflow in MCG-6-30-1 is about 5000 km s _1 using Monte Carlo simulations of the 
K a line propagation in the constant-velocity wind. They present the results of fitting the 
wind model to the XMM (EPIC) data of Wilms et al. (2001) obtained during the observation 
of MCG-6-30-15. These observations were focused on the broad Fe Ka line at ~ 6.4 keV. 
The continuum has been fitted with a power law of index 1.8 and the residuals have been 
fitted using the wind model. The best-fit model parameters they found are E ph = 6.5 keV, 
kT = 0.1 keV, r T = 1.7, (5 = v/c = 0.02 where E ph = hv ph is energy of the primary K a 
photons, and kT is the temperature of the wind respectively. Downscattering of the primary 
line photons in the outflow leads to the formation of red-skewed lines. This effect is seen 
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even when the outflow optical depth r T is about 1-2 and (3>v/c~ 0.02. 

Turolla et al. (1996) were the first to analyze the radiative transfer in expanding (di- 
verging) flow using a Fokker-Planck equation derived by Blandford and Payne (1981). A 
thorough analysis and review of the diffusion theory of photon propagation in an optically 
thick fluid in a bulk outflow has been provided by Titarchuk, Kazanas & Becker (2003), 
hereafter TKB03. They show that the iron line is produced in an effectively optically thick 
medium. Its red wing is the result of multiple scattering, each scattering producing a first 
order v/c— redshift (see Fig.l for a simple explanation of the outflow redshift effect). This 
process produces a red wing to the line without any particularly fine-tuned geometric ar- 
rangement. The TKB03 results are obtained using the Fokker-Planck equation [Blandford & 
Payne (1981)] which is valid for (3 = v/c <C 1 and r e // > 1. It is worth noting that recently, 
using this Fokker-Planck technique, Titarchuk & Shrader (2004) presented a detailed study 
of the outflow downscattering effect in the continuum. The Monte Carlo (MC) method does 
not have these limitations in terms of (3 and r e //. Using a Monte Carlo method Laurent & 
Titarchuk (2004), hereafter LT04, have made extensive radiative transfer simulations in a 
bulk outflow from a compact object (black hole or neutron star). They find that electron 
scattering of the photon within the expanding flow leads to a decrease of its energy which is 
of first order in v/c. LT04 demonstrate that the emergent line iron profile is closely related 
to the time distribution of photons diffusing through the flow (the light curve) and exhibits 
a broad redshifted feature. LT04 find so called negative time lags, related to the time depen- 
dence of the photon energy losses from propagating out of the flow. Thus one can conclude 
that redshifted lines are intrinsic properties of the powerful outflows that are observed in 
many compact objects. 

In this paper we argue that the Fe K a line is naturally produced in an external outflow 
illuminated by the hard X-ray radiation coming from the central source (BH), and the red- 
skewed features are the result of scattering of iron line photons in the outflow. We offer a 
model in which the hard radiation of the central object illuminates and heats the outflow 
region originated in the outskirts of the disk (well outside the innermost part of the accretion 
disk near the BH). We present details of self-consistent calculations of atomic features and 
temperature structure within the ouflow. The result of these calculations for (3 = gives us 
the temperature and ionization structure within the disk. 

In §2 we present a description of the outflow illumination model and provide the main 
idea of the photon frequency shift due to electron scattering in the diverging flow (outflow) 
and converging flow (inflow). We show the results of calculations of temperature and ioniza- 
tion balance in the outflow and evaluation of the equivalent widths of Fe Ka as a function of 
the ionization parameter for various incident X-ray spectral distributions and outflow bulk 
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velocity f3, in §3. Discussion and conclusions follow in §4. 



2. Outflow Illumination by Black Hole X-ray Radiation: Effect of the K a — 
Photon Outflow Redshift and Suppression of the Blue Wing of K a — line 

Our basic scenario is illustrated in Figure 1. The wind originates at a distance Ri nner 
from the central black hole, and is of a density such as to give a Thomson scattering optical 
depth r T close to unity far from the black hole. The optical depth in the Fe K continuum 
is about 1-3 times higher than that due to electron scattering (assuming a solar abundance 
of Fe, and depending on charge state, see Kallman et al. 1994), and so Fe Ka formed by 
inner shell ionization of Fe ions in the outflow by the continuum from the central black 
hole only comes from a smaller inner region. So long as R (rF e K ~ 1) < R(tt ~ 1), the Fe 
Ka equivalent width is insensitive to the Fe abundance assumed. A lower Fe abundance 
increases R (rp e K ~ 1), an d consequently Fe Ka forms over a larger volume. This figure also 
illustrates the redshift of photons in the diverging flow. A photon emitted near the inner 
boundary and subsequently scattered by an electron moving with velocity Vi, is received by 
an electron moving with velocity v 2 as shown with frequency u 2 = + (vi — v 2 ) -n/c] where 
n is a unit vector along the path of the photon at the scattering point. In a diverging flow 
(vx — v 2 ) • n/c < and photons are successively redshifted, until scattered to an observer 
at infinity. In a converging flow (vi — v 2 ) - n/c > and photons are blueshifted. Any 
photon with energy about 7 keV and higher interacting with outflow plasma is more likely 
to be absorbed by the flow and be reemitted at energies about 6.4-6.6 keV (depending on 
the ionization stage of the flow) instead of being scattered by electrons there. This photo- 
absorption effect is particularly important in the view of the main claim of this work as a 
suppression of the blue wing at the expense of the K a emission at energies about 6.4-6.6 keV. 



3. Temperature and Ionization Balance in the Outflow 

Formulation of the Problem 

The temperature and ionization balance balance in the outflow is determined by seeking 
the temperature at which the outflowing gas attains photoionization-recombination equilib- 
rium. The gas is heated by Compton scattering and photoionizations by photons from the 
central compact object, and cooled by radiation, ionization and adiabatic expansion losses. 
We take collisional ionization and recombination rates from Mazzotta et al. (1998), photoion- 
ization cross sections from Verner et al. (1996) and radiative power losses from Summers & 
McWhirter (1979). For the Fe K photoionizations we modify the Verner et al. (1996) cross 
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sections by a multiplicative factor in the range 1-3 depending on charge state to bring these 
cross sections into better agreement with the more recent work of Kallman et al. (2004). It 
is evident that the temperature and ionization structure in the disk due to illumination of 
the hard radiation can be obtained as a result of this solution for f3 = v/c = 0. 

Results of Calculations 

For each inner wind radius, we define a density which gives an electron scattering opac- 
ity in the wind tt, thus density n oc T T /R inner , where Ri nner is the inner radius of the wind. 
In outer portions of the wind, the density varies as n oc l/R 2 , so the ionization parameter, 
L/nR 2 , remains approximately constant. The Compton and photoelectric heating rates, 
approximately proportional to n/R 2 oc l/R 4 , fall off slightly faster than the cooling rate 
(proportional to n/R oc l/R? for adiabatic expansion, proportional to n 2 oc 1/R A for radia- 
tive losses) in a realistic wind. The temperature we calculate at the inner edge of the wind 
is the highest temperature in the flow in cases where adiabatic expansion dominates the 
cooling. Where radiative cooling is more important, the heating and cooling are more evenly 
balanced, and the temperature we calculate is likely to be approximately correct through- 
out the flow. Detailed radiation transfer modelling accounting for successive absorption and 
scattering of the incident hard spectrum to improve upon this estimate is beyond the scope 
of this work. 

In Figure 2 we plot the run of plasma temperature against initial wind radius for a variety 
of models. The top panel shows temperature (solid line, to be read on the left hand y-axis) 
and average Fe charge state (dashed line, to read read on the right hand y-axis) for an outflow 
with (3 = v/c = 0.1 and for incident spectra I{y) oc v~ a exp(— hv/2kT) with a = 0.25, 0.5, 
0.75 and kT = 50 keV. The inverse of the ionization parameter nR 2 /L ~ T^R inner /L m on 
this plot, where L 40 is the source luminosity in 10 40 ergs s" 1 . For T^Ri nner / L m < 10 12 ' 8 cm, 
Fe is in the He-like charge state or higher, and adiabatic expansion dominates the cooling. 
As L-shell ions Fe XVII- XXIV start to form, radiative cooling begins to dominate and the 
temperature drops more steeply, with regions of thermal instability between Fe charge states 
16+ and 13+. Harder spectra (i.e. lower a) with more photons at high energies have slightly 
lower photoionization rates for charge states Fe 23+ and below (for the same luminosity), and 
so under these conditions the outflowing plasma recombines and cools slightly faster. The 
lower panel shows similar plots, but now with a = 0.5 and (3 = v/c = 0.02, 0.05, 0.1, and 
0.25. Faster outflows give more adiabatic cooling and so cool and recombine at lower radii 
than do the lower outflows. 

Figure 3 shows the equivalent widths of Fe Ka (dashed lines, to be read on the right 
hand axis) for the same sets of parameters as for Figure 2. The solid lines to be read on the 
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left hand axis give the total opacity in the outflow at 4 keV, in the form exp (—opacity). At 
small radii the equivalent width is essentially zero because Fe is in charge state 24 or higher, 
with no L shell electrons available to fill a K shell vacancy. As Fe recombines at larger radii, 
the equivalent width increases. The predicted equivalent width at large inner radii of about 5 
keV is significantly larger than that observed of about 1 keV. However the total wind opacity 
at these radii is such that most of the Fe K line would be absorbed further out. Only for 
Rinner ~ 10 13 L 40 and consequently for the wind temperatures about a few times 10 6 K (see 
Figs. 2-3) to where the equivalent width is approximately 1 keV is the opacity sufficiently 
small to allow the line to be observed. We also note that this opacity is likely to have the 
effect of suppressing any Fe L shell features between 1 and 2 keV photon energy. 



4. Discussion and Conclusions 

We find the self-consistent temperature and ionization structure of the wind shell as a 
function of the parameter (radius/luminosity, T^R inner / 'L 40 , where Ri nner is the inner outflow 
radius and L 40 is X-ray luminosity in units of 10 40 erg s -1 ) for a given shell Thomson optical 
depth tt- It is evident that the ionization parameter L^/nR 2 is constant through the shell 
(i.e. for R > Ri nner ) if the velocity of the wind is constant through the flow. The constraints 
on the Fe ionization balance and the opacity in the wind dictate a possible range of interest 
of TTRi nner / X40 values around 10 13 cm. Thus our solution allows to determine the size of 
the shell base for a given luminosity for which the red-skewed line is observed. For a typical 
low/hard-state luminosities of a few times 10 36-37 erg s _1 when the red-skewed line is really 
observed in galactic black holes the inner radius of the shell is a few times 10 9 ~ 10 cm for 
the optical depth of the wind of order unity. It means that the inner radius of the wind 
shell is about 10 3-4 of Schwarzschild radii. One can obtain the similar size of the wind 
shell in Schwarzschild radii units for extragalactic sources because the luminosity and the 
Schwarzschild radius is linearly scaled with BH mass. 

It is important to note that if Thomson optical depth of outflow is of order of unity then 
the mass outflow rate M out is of order of the Eddigton mass accretion rate M E dd and higher 
(see e.g. formula 4 in King & Pounds 2003, hereafter KP03). King (2003) and KP03 present 
strong arguments that powerful mass outflows from Eddington-limited accreting compact 
objects appear to be a very widespread phenomenon. They further argue that they may 
provide the soft excess observed in quasars and ULXs, and imply that such objects have 
a major effect on their surroundings. Recent XMM-Newton observations of bright quasars 
(Pounds et al., 2003a,b; Reeves et al., 2003) give strong evidence for powerful outflows 
from the nucleus with mass rates M out ~ M Q yr -1 ~ MEdd and velocity v ~ 0.1c (i.e. 
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(3 = 0.1) in the form of blueshifted X-ray absorption lines. These outflows closely resemble 
those resently inferred in a set of ultraluminous X-ray sources with extremely soft spectral 
components (Mukai et al. 2003; Fabbiano et al., 2003). Furthermore, if the outflow optical 
depth is really about one in the low/hard state (when the X-ray luminosity is much less the 
Eddington luminosity) one can conclude that in the low/hard state the disk mass accretion 
rate is only very small fraction of the outflow mass rate. In fact, this ratio of the accretion 
and outflow rates was predicted by Blandford and Begelman (1999) who developed the pure 
hydrodynamical model (so called "ADIOS" model). The main point of the ADIOS model is 
that only a tiny fraction of the gas supplied actually falls on to black hole. This is precisely 
what we obtain using the diverging flow model. 

We conclude that the range of the parameter T^R inner / 'L 40 (proportional to the inverse 
of the so called "ionization parameter" used in the literature) is about 10 13 when the observed 
K a lines of the ~ 1 keV equivalent widths are presumably produced in the wind. The wind 
temperature (in such a case) is about a few times 10 6 K. Thus our study shows that the 
strong iron line with its red-skewed feature can be generated in the relatively cold extended 
region far away from the source of the illuminating photons (of order 10 3 ~ 4 of Schwarzschild 
radii). 

From the other hand one can argue that the "standard" cold disk plus hot corona 
explanation for the broad lines is "very robust (it needs only the cold disk and a hot corona)" 
and one does not need any other explanation for this effect. In fact, the so called "cold 
disk" can survive under the hot corona only if the corona is situated at least 10 3 — 10 4 /rx 
Schwarzschild radii where tt is the Thomson optical depth of the disk. If corona is much 
closer the disk as a target illuminated by the hard radiation of the corona is very hot and 
finally evaporated. It is evident from Figs. 2-3, that temperature goes to 10 9 K and and 
K a equivalent width goes to zero as T T R inner / X 40 <C 10 13 cm for (3 < 0.05. Other points 
regarding the "robustness" of the standard interpretation are: even the line is produced in 
the disk how one can see K a through the corona if the corona optical depth r cor is of order 
of one and the corona electron temperature is around 60 keV? The directed radiation of the 
K a emission is attenuated exponentially as exp(— r cor //i) (at least) < 1/3 where fi is the 
cosine of the inclination angle. But the scattered component is completely smeared out - the 
relative energy change of the line photon at any scattering is < AE/E >= AkT/m e c 2 = 0.5 
(see e.g. Sunyaev & Titarchuk 1980). But even if the disk line photons manage to come to 
the observer one would see them with very prominent blue wing, which is a result of first 
order v/c effect, as the gravitational redshift is an effect of second order with respect to v/c, 
i.e Rs/R ~ ( v / c ) 2 where R$ is Schwarzschild radius. The demonstrated application of our 
outflow model to data points out a potential powerful spectral diagnostic for probes of the 
outflow-central object connection in Galactic and extragalactic BH sources. 
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Fig. 1. — On the left side: Schematic diagram showing wind geometry. The outflow (wind) 
originates at the inner radius. The radii where the optical depth in the Fe K continuum 
and the electron scattering optical depth are of order unity. A photon emitted near the 
inner boundary and subsequently scattered by an electron moving with velocity vi, is 
received by an electron moving with velocity V2 as shown. The change in frequency is 
^2 = ^i (1 + (vi — V2) • n/c) where n is a unit vector along the path of the photon scattered 
at (b). In a diverging flow (vi — v 2 ) ■ n/c < and photons are successively redshifted, until 
scattered to an observer at infinity. The color of photon path indicates the frequency shift 
in the rest frame of the receiver (electron or the Earth observer). On the right side: In a 
converging flow (vi — v 2 ) • n/c > and photons are blueshifted. 
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Fig. 2. — Plot of temperature at Rinner against Rinner for an accretion wind (solid lines). 
The inner radius Ri nne r scales with the the central source luminosity and the optical depth 
of the wind shell tt- For any r T ^ 1 the temperature is a function of TiR inner / L 40 . The top 
panel shows results for expansion speed /3 = v/c = 0.1 and spectral indices a = 0.25, 0.5, 
and 0.75. Bottom panel shows results for a = 0.5 and (3 = v/c = 0.02, 0.05, 0.1, and 0.25. 
The dashed lines show the average Fe charge state (to be read on the right hand axis). 
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Fig. 3. — Plots of the equivalent width of Fe Ka (dash lines) produced by inner shell 
ionization of Fe in an accretion wind by continuum from the central source, for the same 
parameters as is Figure 2. At small radii the equivalent width is essentially zero because Fe 
is in charge state 24 or higher, with no L shell electrons available to fill a K shell vacancy. As 
Fe recombines at larger radii, the equivalent width increases. The solid lines to be read on 
the left hand axis give the total opacity in the outflow at 4 keV, in the form exp (—opacity) . 
The predicted equivalent width at large radii of about 5 keV is significantly larger than that 
observed of about 1 keV. However the total wind opacity at these radii is such that most 
of the Fe K line would be absorbed further out. Only for T rm TRi nn er/ £40 out to where the 

equivalent width is approximately 1 keV is the opacity sufficiently small to allow the line to 
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